Measuring acute effects of subanesthetic ketamine on cerebrovascular hemodynamics in humans using TD-fNIRS

Quantifying neural activity in natural conditions (i.e. conditions comparable to the standard clinical patient experience) during the administration of psychedelics may further our scientific understanding of the effects and mechanisms of action. This data may facilitate the discovery of novel biomarkers enabling more personalized treatments and improved patient outcomes. In this single-blind, placebo-controlled study with a non-randomized design, we use time-domain functional near-infrared spectroscopy (TD-fNIRS) to measure acute brain dynamics after intramuscular subanesthetic ketamine (0.75 mg/kg) and placebo (saline) administration in healthy participants (n = 15, 8 females, 7 males, age 32.4 ± 7.5 years) in a clinical setting. We found that the ketamine administration caused an altered state of consciousness and changes in systemic physiology (e.g. increase in pulse rate and electrodermal activity). Furthermore, ketamine led to a brain-wide reduction in the fractional amplitude of low frequency fluctuations, and a decrease in the global brain connectivity of the prefrontal region. Lastly, we provide preliminary evidence that a combination of neural and physiological metrics may serve as predictors of subjective mystical experiences and reductions in depressive symptomatology. Overall, our study demonstrated the successful application of fNIRS neuroimaging to study the physiological effects of the psychoactive substance ketamine in humans, and can be regarded as an important step toward larger scale clinical fNIRS studies that can quantify the impact of psychedelics on the brain in standard clinical settings.

www.nature.com/scientificreports/ measurements, including a photoplethysmography sensor (PPG), an electrodermal activity sensor (EDA), and a nasal cannula to measure end-tidal CO 2 (EtCO 2 ) (Fig. 1d). The amount of administered ketamine was determined by the participants' weight with a target of 0.75 mg/kg, not to exceed a total dose of 60 mg (reached by two participants) (Fig. 1e).
Ketamine induces altered states of consciousness and reduces depressive symptomatology. Participants completed surveys at different time points in the course of the entire study ( Fig. 1a; "Methods"; Supplementary Fig. 1), which served to evaluate the subjective (phenomenological) effects of the ketamine vs. saline injection. All participants had higher total scores on the Revised Mystical Experience Questionnaire (RMEQ; paired t-test, p < 10 −8 ), the 5-Dimensional Altered State of Consciousness Rating Scale (5D-ASC; paired t-test, p < 10 −6 ), and the Clinician Administered Dissociative States Scales (CADSS; paired t-test, p < 10 −5 ) after ketamine vs. saline. Most participants also had a higher total score on the Brief Psychiatric Rating Scale (BPRS; paired t-test, p < 10 −3 ). In addition, participants completed the Quick Inventory of Depressive Symptomatology Ketamine affects the cardiovascular system, absolute cerebrovascular hemoglobin oxygenation and concentration but preserves respiratory parameters. Often regarded as artifacts, fNIRS data are rich in physiological information from which multiple physiological features can be extracted 48,49 . Although ketamine is thought to preserve baseline respiratory parameters 50 , it is sympathomimetic and has known effects on cardiac signals such as heart rate (HR) and heart rate variability (HRV) 51 . In fact, cardiac measures alone have been suggested as potential biomarkers of treatment outcomes in patients with major depressive disorder (MDD) 51,52 . Consequently, we tested our ability to retrieve several physiological variables from our data and further demonstrate the effects of the active drug as compared to the placebo control. First, we quantified changes in cardiac signals following ketamine and saline injection. Because this is the first study in which the Kernel Flow1 TD-fNIRS headset was used to derive these signals, we used an auxiliary PPG sensor (attached to the finger) to validate the procedure. Both the PPG sensor and the Flow1 TD-fNIRS system measure pulse rate (PR), which is strongly related to HR 48 . A periodic pulse-like signal (representing the blood volume pulsation) is very apparent when visualizing total photon counts from the 850 nm wavelength of the TD-fNIRS system averaged over well-coupled prefrontal channels, which is aligned with the blood volume pulsation measured by the PPG sensor (Fig. 2a). We found that there was a strong agreement between the instantaneous pulse rate (PR) computed from the PPG and that from the Kernel Flow1 TD-fNIRS system ( Fig. 2b; Pearson correlation r = 0.93 ± 0.02; mean ± standard error of the mean (SEM) over all sessions) ("Methods"). We also compared the average PR and pulse rate variation (PRV) extracted from the PPG and the Flow1 system and found that PPG-derived and Flow1-derived metrics were strongly correlated (Fig. 2c). We proceeded to characterize the changes in PR and PRV following ketamine and saline administration. We found a significant increase in PR during the ketamine dosing session (82%) (Fig. 2d) compared to the saline session as well as a significant decrease in PRV (91%) (Fig. 2e)-both findings were confirmed with the external PPG sensor as well and are consistent with those described in prior literature 51 .
Next, we asked whether there were significant differences between ketamine and saline dosing sessions with respect to cerebrovascular hemodynamics and oxygenation. Therefore, we computed the absolute concentrations of oxy-and deoxy-hemoglobin concentrations commonly referred to as HbO and HbR, respectively ("Methods"). To verify that these signals exhibited expected physiological patterns, we studied the breathing challenge performed during the baseline run prior to each dosing session. Here, a brief hypocapnic state was induced at minute 8 as participants were instructed to breathe rapidly. Indeed, absolute HbO and HbR time-series showed the expected change during the hypocapnic state ( Supplementary Fig. 2a). To quantify the change in hemoglobin concentrations as a result of saline or ketamine administration, the session-wise median value of HbO and HbR time-series during the resting-state baseline run were subtracted from those during the dosing sessions. We found a significant increase in HbO (paired t-test, p = 0.032) and a significant decrease in HbR (paired t-test, p = 0.021) during the ketamine administration as compared with saline sessions (Supplementary Fig. 2b).
Lastly, we explored other physiological measures that were captured in parallel. Using the capnography timeseries recorded from the Nonin capnograph ( Supplementary Fig. 3a), we computed the breathing rate during saline and ketamine dosing sessions, and found no significant differences between the two sessions (p > 0.1; Supplementary Fig. 3b; "Methods"). Similarly, no significant differences were observed when comparing the average arterial oxygenation (SpO 2 ) and EtCO 2 values from ketamine sessions to those from saline sessions ( Supplementary Fig. 3b), as expected 50 . We did, however, observe an increase in the electrodermal activity (EDA) during ketamine sessions as compared to saline sessions (Supplementary Fig. 3c; saline: 0.19 ± 0.02 µS, ketamine: 1.51 ± 0.52 µS; mean ± SEM). In addition to the subjective measures of experience following ketamine administration, significant changes in PR, PRV, and EDA provide further corroboration for the effectiveness of the ketamine doses administered in the study.
Ketamine reduces the whole brain fractional amplitude of low-frequency fluctuations. The fractional amplitude of low-frequency fluctuations (fALFF) in the fMRI-derived blood-oxygen-level-dependent (BOLD) signal has been shown to correlate with cerebral metabolism and is a proxy for local neuronal activity 47 . Abnormal fALFF has been implicated in a variety of conditions, including unipolar and bipolar depression [53][54][55] . Of note, fALFF within the default mode network (specifically hypoactivity of slow 5 fALFF) is implicated in depression and alleviated by treatment. An acute reduction in fALFF following ketamine administration has been reported previously using fMRI 56 . To investigate the changes in fALFF in our study, we first converted the distributions of the times of flight of photons (DTOFs) to the relative concentrations of HbO and HbR ("Methods"). We then calculated the whole-brain time-varying fALFF for both chromophores using a sliding window approach (duration: 5 min, stride: 10 s; "Methods"). We found that, at the group level, ketamine dosing sessions showed a marked decrease in fALFF (Fig. 3b), in particular in the later part of the session. Thus, we computed whole-brain fALFF for the first half and second half of the dosing sessions separately. To account for any baseline differences, we used the pre-breathing-challenge portion of the preceding resting-state baseline run to compute a baseline fALFF ("Methods"). This time interval, lasting for 400 s, was considered to be a resting session (as it was prior to the breathing challenge) and was considered as a "baseline" segment (i.e. no stimulus-or taskevoked physiological activity) here and throughout the analyses when baselining was warranted. Baseline-normalized whole-brain fALFF was significantly lower in the second half of the ketamine dosing sessions compared to the second half of the saline dosing sessions, for both HbO (paired t-test, p = 5.97 × 10 −3 ) and HbR (paired t-test, p = 0.025) (Fig. 3c) www.nature.com/scientificreports/ on whole-brain fALFF reaches its maximum around the same time as the PR peaks (at the 500 s mark, on average; Fig. 3a); however, the reduction in fALFF appears to be more sustained than PR, with decreased fALFF persisting for at least 30 min after the injection, at a time when the PR has returned to baseline. This suggests that changes in fALFF and the cardiovascular system are independently affected by the ketamine injection.
Ketamine acutely reduces global brain connectivity of the prefrontal cortex. The TD-fNIRS system used in this study provides whole head coverage and is thus well-suited to study patterns of cortical connectivity and their acute changes following ketamine administration. Existing fMRI studies focusing on ketamine-induced acute connectivity changes in healthy volunteers found an increase in brain connectivity 19 . However, this finding is not without controversy. For example, using the measure of global brain connectivity (GBC), some studies report the detection of a global increase in functional connectivity (in particular in prefrontal cortices) [57][58][59][60] , but an independent study from another group failed to replicate this finding 61 . Studies using seedbased connectivity report some up-and some down-regulated connections: for example the connection between the posterior cingulate cortex and medial prefrontal cortex is typically found to be downregulated 62,63 ; decreased connectivity has also been reported in the visual network 64 ; increases and decreases have been reported between the dorsolateral prefrontal cortex and the rest of the brain 65 . In light of the extant literature, our expectations were therefore rather mixed.  www.nature.com/scientificreports/ To calculate functional connectivity, TD-fNIRS data (i.e. HbO and HbR time-series) underwent further preprocessing ("Methods"). The group-averaged "dense", channel-space connectivity matrix ( Fig. 4a) exhibited hallmark features of whole-brain functional connectivity such as community structure and homotopic connections. We computed connectivity matrices separately for each session of each visit: baseline and dosing sessions, during the saline and the ketamine visit. For dosing sessions, we focused on a 400 s segment (≈ 6.5 min) at the end of the session, leveraging the observation that fALFF showed the largest differential between ketamine and saline in the later part of the session. The duration of this segment was chosen such that it was the same length as the baseline segment to avoid introducing potential confounds. We were primarily interested in comparing the ketamine and saline dosing sessions. As with the fALFF analysis, we corrected for any differences in baseline connectivity between the saline and the ketamine dosing sessions by using the baseline segment (400 s from the immediately preceding baseline session; "Methods")-such baseline differences, whether accidental or related to participants' expectations, may obscure results or, alternatively, mislead their interpretation.
We started by assessing GBC, a graph-theory inspired metric (a.k.a. weighted degree centrality), which summarizes the connectivity strength of each node of the connectivity matrix 66,67 . Because the prefrontal area is Figure 3. Changes in whole-brain fALFF following ketamine administration. (a, b) Time-varying PR and fALFF were normalized by the initial value at a given session to obtain the percentage of change throughout the session. Shown are changes in PR (a) and whole-brain fALFF for HbO and HbR (b) during saline (gray) and ketamine (green) dosing sessions (mean ± SEM). (c) During each saline or ketamine session, whole-brain fALFF (i.e. averaged across the brain) was computed for the first half and second half of the sessions independently and was baseline normalized. Only during the second half of the sessions, fALFF values were significantly lower for the ketamine session compared to the saline session, for both HbO and HbR (paired t-test, p = 5.97 × 10 −3 and p = 0.025 respectively). www.nature.com/scientificreports/ where GBC has been found to be most affected in previous studies 57-60 , we focused our analysis on this region (i.e. quantifying the connectivity of prefrontal channels to the rest of the brain)("Methods"). As explained above, we baselined GBC for each dosing session using GBC computed during the pre-challenge portion of the restingstate baseline session collected in the same visit. In this setting, we found that ketamine acutely reduces prefrontal GBC compared to saline ( Fig. 4b): though no channel or cluster of channels reached significance on its own, the distribution of prefrontal GBC at the group level was significantly shifted to the left (indicating reduction, Fig. 4c). When the participants' prefrontal GBC was collapsed (by taking the median value of each individual's distribution over prefrontal channels), indeed, we found a significant reduction in GBC for HbO (p = 8.8 × 10 −3 ) and a trend towards reduction for HbR (p = 0.05). We further investigated whether any connections were downregulated or upregulated by ketamine, compared to saline: we did not find any connections that were significantly more modulated by ketamine than by saline injection ("Methods").
Physiological and neural features and their relationship with subjective mystical experiences and depressive symptoms. It is argued that one of the predictors of outcome measures in psychedelic therapy is the intensity of mystical experiences reported by the participants 68 . As such, exploring parameters capable of such predictive power can pave the way for treatment optimization. Could a combination of neural and physiological features predict the degree of participants' mystical experiences? To answer this, we regressed participants' self-reported RMEQ total scores against time-varying fALFF and PR features (Fig. 5a, Leave One Figure 5. Changes in the magnitude of whole-brain fALFF following ketamine administration. (a) Timevarying PR and fALFF (normalized by the initial value at a given session) exhibited a separation between ketamine (green cluster) and saline (gray cluster) sessions (note the ketamine cluster in the bottom right corner indicating an increase in the PR and decrease in fALFF). (b) By using a linear regression model with time-varying PR and fALFF (both HbO and HbR) as input features, we were able to predict the total RMEQ score. Shown are the model predictions on the validation set during each fold (mean ± SEM of the timevarying predictions within each session versus the reported RMEQ total score for the entire session). (c) We evaluated the trained model on the time-varying features to obtain time-varying RMEQ score predictions (on the validation set). Shown are the normalized prediction scores sorted by the peak time across participants. It appears that different participants may reach their peak of mystical experiences at different time points after ketamine administration. (d) GEE model revealed that the interaction term between session type (ketamine vs. saline) and whole-brain fALFF (whether it was high or low throughout the session) was a significant factor in the change in QIDS scores (HbO: p = 3.53 × 10 −6 , HbR: p = 4.76 × 10 −7 ). Shown are the change in QIDS scores after the ketamine sessions when we performed a median split on the participants' fALFF during their ketamine session. www.nature.com/scientificreports/ Subject Out cross-validation technique; "Methods"). The predicted scores on the held-out sessions showed a mean squared error of 1.77 ± 0.44 (mean ± SEM, see Fig. 5b), and the model r-squared on the held-out set was 0.24. Furthermore, we ensured that the performance of the model was above chance level by re-computing the results for 1000 random permutations of the RMEQ scores (p = 0.004; "Methods"). We also calculated the partial correlations between each feature and the RMEQ score while using the other features as covariates, to tease apart the contributions of cerebrovascular and systemic physiological features. We found PR, HbO fALFF, and HbR fALFF to have similar full correlations (0.57, − 0.47, − 0.43, respectively) as well as partial correlations (0.49, − 0.36, − 0.32, respectively; p < 10 −5 for all), suggesting that cerebrovascular and systemic physiological signals contain complementary and non-redundant information. When the full model was evaluated at different time points within the ketamine dosing session to obtain time-varying predicted RMEQ scores, it appeared that the peak prediction occurred at different times for each participant (Fig. 5c).
Lastly, we asked whether the change in QIDS scores were related to fALFF during the session. To this end, we modeled the change in QIDS scores as a function of fALFF during dosing sessions (in particular whether session fALFF was high or low, as defined by a median-split), session type (saline/ketamine), and the interaction between the two terms using Generalized Estimating Equations (GEEs). We found that the change in depressive symptoms not only depended on ketamine vs. saline but also on fALFF such that participants with higher fALFF showed a larger reduction in depressive symptoms after ketamine ( Fig. 5d; interaction term p < 0.001 for both HbO and HbR). Whether the combined changes in PR and fALFF can provide therapeutically relevant biomarkers remains to be explored in future large scale studies in patients.

Discussion
In the current study, we report on behavioral, systemic physiological, and cerebrovascular measurements from participants who underwent neuroimaging using TD-fNIRS while experiencing altered states of consciousness as a result of subanesthetic ketamine injection in a clinical environment. We compared the acute effects of intramuscular bolus ketamine administration against a placebo (saline administration) condition. We demonstrate reliable physiological measurements of PR and PRV extracted from TD-fNIRS recordings that match those obtained from commercial PPG sensors, thus rendering the use of external sensors to measure cardiac activity unnecessary in future experiments. We also find robust differences between ketamine and saline conditions: (1) physiological measures: increased PR, decreased PRV, increased absolute HbO and decreased HbR, and elevated EDA (measured by an external sensor) during ketamine dosing, and (2) neural measures: decrease in fALFF and in the global brain functional connectivity of the prefrontal cortex following ketamine administration. These results are generally consistent with findings reported by previous studies (apart from the global functional connectivity result) 51,56,69,70 .
By using a combination of systemic physiological (i.e. PR) and neural (i.e. fALFF) features, we were able to predict participants' reported RMEQ scores. Further, we found that when the trained model was evaluated throughout the ketamine session, the predicted peak of mystical experience occurred at different times for different individuals. The predicted time-series approximate the expected pharmacokinetics after such an intramuscular bolus and may reflect the time course of the psychoactive effects of ketamine in each participant; to corroborate this, future studies could include continuous measurements of plasma concentration of the psychoactive substance as a proxy for continuous subjective experience 71 . Subjective experiences, particularly those along the dimensions of anxiety and mystical experience, are of particular interest as they have been suggested to affect treatment outcome measures when ketamine is used as an antidepressant 72,73 . Thus, exploring the relationship between the different dimensions of subjective experience and neural data is an important question for future research.
Although we observed a decrease in whole-brain fALFF following ketamine administration, this effect may be region-or network-specific. It has been suggested that fALFF within the default mode network (DMN) is of particular functional importance 54 . Indeed, individuals with MDD show increased fALFF in the left dorsomedial prefrontal cortex, and a ketamine-induced decrease in fALFF suggests a plausible mechanism by which ketamine alleviates depression symptoms 27,74,75 . Furthermore, we found that fALFF during the ketamine session was a significant factor in the change in depressive symptoms following ketamine administration. It must be borne in mind, however, that while fALFF may be a useful biomarker in predicting treatment outcomes, our study involved only healthy volunteers, and caution must be exercised when extending this measure to the clinical population 27 .
Our result of acutely decreased global functional connectivity in the prefrontal cortex after ketamine administration appears to be at odds with extant literature 59,60 (where a slightly lower dose was used as in our study [0.5 mg/kg vs. 0.75 mg/kg]) though there was a recent report of a failure to replicate past positive findings of a GBC increase in healthy volunteers 61 . Despite the focus on prefrontal GBC in the current manuscript and prior literature, future large scale studies are needed to elucidate whether this effect is limited to the prefrontal region or if it can be observed in other regions as well. In addition, our failure to identify any specific connections with significantly modulated connectivity is at odds with a limited, but certainly mixed, body of reports of acute functional connectivity alterations following ketamine injection. Current studies, this included, have been small and large-scale studies are needed to draw definitive conclusions. We also acknowledge that there are inherent limitations to the whole-head functional connectivity estimates derived with the TD-fNIRS prototype used in this study: for example, not all participants had perfect coupling over the whole head, and the registration of data across participants can be further improved, which may increase statistical power at the group level. Future work is required to ensure improved measurement coverage over the head, and better interparticipant registration capabilities. In this work, reported measures are averaged over large anatomical areas (the whole head or the entire prefrontal cortex), however more sophisticated anatomical modeling and registration methods may allow for more detailed investigations of brain measures. In addition, there are differences in the www.nature.com/scientificreports/ neuronal-activity-related cerebrovascular changes measured with fMRI and fNIRS. fNIRS is more sensitive to changes in the venous blood fraction and capillaries, compared to the BOLD fMRI signal, and both fMRI and fNIRS measures are a proxy for neural changes 76 . Consequently, systemic physiological activity can have a significant impact on fNIRS signals 77,78 , although our measurement and signal-processing approach was optimized to reduce this impact. The current study, due to its small sample size and limitations of only measuring cortical connectivity, was not designed to adjudicate conflicts in the existing literature. It was intended as a proof-of-principle for using a novel TD-fNIRS system in clinical settings. The simplicity of using TD-fNIRS in the clinic, and its ability to measure cortical function (fALFF, functional connectivity) in longitudinal and large-scale studies, may unlock insights on functional brain reorganization, acutely and chronically, following pharmacological interventions. These new insights may be particularly valuable for tailoring KAP or extended to tailoring PAP, more generally.
This study has several further limitations. Our study population was not only small, it was also a sample of healthy volunteers rather than a patient population. If the intent is to understand the mechanism of action of RAADs in clinically depressed patients, a MDD patient sample would be most appropriate. Thus, the results that we report in this work may not apply directly to ketamine treatment of depressed patients. Another limitation is that, while we measured subjective depressive symptomatology one week after ketamine dosing, we did not include brain measurements post-dosing which is when many studies have noted neural changes and psychological effects are measurable (see 19 for a review of acute vs. delayed effects of ketamine administration). Our design was also single-blind (instead of double-blind) and used a fixed predefined order of sessions (saline first and ketamine second). Even with this order, because the drug induces such potent subjective experiences, most but not all participants were effectively unblinded after participating in the first (saline) dosing session (Supplementary Fig. 1). To avoid completely unblinding the participant after their first visit, we did not counterbalance the session order and, therefore, cannot take into account potential order effects. Such issues with blinding are prevalent amongst studies utilizing full doses of psychedelics, and it appears ketamine is no different. Of course, all these shortcomings are not unheard of in a pre-clinical, proof-of-principle study.
Whole head TD-fNIRS neuroimaging offers an ecologically valid methodology to measure the acute and delayed physiological and neural changes following pharmacological interventions, such as ketamine administration. This scalable technology that captures high-quality cerebrovascular signals and can be deployed widely in clinics may be the solution to the current lack of consistency in neuroimaging investigations of the effects of ketamine administration. But why prefer cerebrovascular measurements to electrophysiological measurements, such as those that can be captured with EEG? While EEG suffers from inherent limitations such as its sensitivity to artifacts (muscle, motion, ambient electromagnetic fields) and lack of spatial specificity, hemodynamic estimates of neural activity are themselves plagued with contamination from systemic physiology 77,78 and, in the case of fNIRS, limited depth sensitivity. Future studies may benefit from combining multiple imaging modalities (e.g. EEG with TD-fNIRS) to provide a more complete understanding of the impact of psychedelics on the brain in a clinical setting.

Methods
Participants and screening procedures. Fifteen healthy participants (8 females and 7 males, age (mean ± SD): 32.4 ± 7.5 years, range 24-48 years) completed four study visits. Participants gave written informed consent before beginning the study in accordance with the ethical review of the Advarra IRB (#Pro00059548), which approved this study, and the Declaration of Helsinki. Pre-screening eligibility determination included the following inclusion criteria: adults between the ages of 21-50; in good general health (as ascertained by medical and psychiatric history and physical examination including vital sign measurements); prior experience (at least once but not more than 20 times) with a classical psychedelic drug (e.g., LSD, psilocybin, mescaline, etc.) or ketamine (if known, no more than seven previous 60 mg daily doses) as confirmed by the study clinician; body mass index between 18 and 30 kg/m 2 ; fluency in English; and ability to provide informed consent. Exclusion criteria included: known history of severe neurologic, cardiovascular, or respiratory disease; known history of hepatic and renal disease and/or insufficient hepatic and renal function (as evidenced by clinically significant abnormal lab results); elevated blood pressure; recent or current substance use disorder (except nicotine); history or current severe psychiatric disorder or condition (including psychosis, schizophrenia, or bi-polar disorder, or study clinician-assessed severe mental illness, such as severe untreated depression); suicidal ideation or behavior per self-report and as assessed by the study clinician at every study visit; recently or currently taking MAOI, stimulant, mood stabilizer, benzodiazepine, or other central nervous system depressant; known allergies to ketamine; or pregnancy or nursing.
Participants received monetary compensation. Once participants consented to participate in the study, they underwent a physical exam, medical history, psychiatric history, and assessment for suicide ideation by the study clinician. Vitals, height, and weight were also measured, and demographic information was collected. Participants then did a short breathing challenge while wearing the Kernel Flow1 TD-fNIRS headset to simulate what future study visits would feel like. Lab tests (complete metabolic panel, pregnancy test, and drug screen) were completed and reviewed by the study clinician no more than 14 days before the first dosing visit. A total of 15 participants completed all of the study visits. Participants were given the option for an additional meeting with the study clinician to ask more detailed questions about what to expect and how to prepare for dosing visits.

Study design.
The study was a single-blind, placebo-controlled, non-randomized design with participants completing study visits roughly once a week for four weeks. The four study visits were always conducted in the same order: a screening visit, two dosing visits, and a follow-up phone call. Dosing visits were always pla- www.nature.com/scientificreports/ cebo (saline, 0.9% NaCl) first and ketamine second, with the ketamine visit occurring one week (7.1 ± 0.5 days, mean ± SD) after the saline visit (Fig. 1a). Ketamine and saline were administered via bolus intramuscular injection (deltoid muscle). Ketamine dosing was based on participant weight with a target of 0.75 mg/kg, up to the maximum dose of 60 mg. Two participants were administered the maximum dose (Fig. 1e).
Behavioral and experimental paradigms. At each dosing visit, participants underwent physiological and brain monitoring during a baseline run featuring a short breathing challenge and for 30 min after saline or ketamine was injected. Physiological sensors (E4, Empatica, Boston, MA, USA and Lifesense, Nonin, Plymouth, MN, USA) measured respiration rate, SpO 2 , EtCO 2 , PR, and EDA. Brain data was measured with a whole-head TD-fNIRS system 33 (Flow1, Kernel, Culver City, CA, USA) (for a detailed description see "Methods": fNIRS data collection and analyses).
During the baseline run, participants were verbally instructed to breathe normally for the first 8 min. Then they were asked to breathe rapidly and deeply until their EtCO 2 dropped by 10 mmHg (hypocapnic state), which normally took 30 s or less. Lastly they were asked to resume breathing normally for 5 min. The entire run was approximately 14 min.
During the saline and ketamine administration, participants were partially reclined on a couch, supported by pillows. The pilot participant demonstrating the setup in Fig. 1b provided written consent for this photo to be shared in an online open access publication. Participants were given the option to listen to music, which was selected by the clinician and was alyrical as well as melodic (14/15 listened to music during both visits, 1/15 declined music at both visits). Room lights were dimmed. Brain and physiological signals were recorded for 30 min post-dose. Vital signs were monitored continuously until 60 min post-dose, and then checked at 90 min and 120 min for safety. Adverse events were recorded as needed throughout the visit. Adverse events occurred in 1 out of 15 (7%) and 9 out of 15 participants (64%) after saline and ketamine injection respectively and included fatigue, high blood pressure, nausea, dizziness, headache, fogginess, thirst, and low SpO 2 ( Supplementary Fig. 4). No severe adverse events occurred. Participants were allowed to leave with an arranged ride a minimum of 120 min post-dose, when deemed safe by the study clinician. fNIRS data collection and analyses.

Data acquisition
Data acquisition was performed using previously described methods 33 . Briefly, the TD-fNIRS system was used to acquire the distributions of the times of flight of photons (DTOFs) from more than 2000 channels across the head. The system consists of 52 modules recording at an effective imaging rate of 7.14 Hz and uses two lasers with wavelengths of 690 nm and 850 nm. Each module contains a source and six surrounding detectors (at 10 mm from the source) and, therefore, multiple channels (a given combination of source-detector) can be formed at various source-detector separations (SDS). The system provides bilateral coverage of prefrontal, parietal, temporal, and occipital regions on separate plates for each region (referred to as: prefrontal, sensorimotor, auditory, and visual respectively) and each hemisphere (Fig. 1c).
In addition to the DTOFs, the Flow1 system is capable of capturing parallel streams of data including: data from 5 inertial measurement units (IMU) distributed throughout the headset as well as temperature sensor data for both sources and detectors that are acquired at all times. Data from auxiliary sources, if connected, will be collected as well (Fig. 1d). For example, in the current study, data from an external photoplethysmography (PPG) sensor was acquired and later used for physiological signal validation.
For the fifteen participants who completed all four study visits, technical issues occurred during some sessions leading to partial (3 of 30 dosing sessions) or full (1 of 30 dosing sessions) data loss ( Supplementary  Fig. 5a).

Data preprocessing and feature extraction: relative hemoglobin concentrations
We implemented a channel selection procedure using both a histogram shape criteria as well as the presence of physiological signals 33 . The percentage of retained channels refers to the number of usable channels of neural data and can vary based on the contact between the headset and participants' head. This percentage is reported for the participants when considering sex, skin type, or hair color and, for the vast majority of the participants, the percentage of retained channels fall between 40 and 80% ( Supplementary Fig. 5b).
Next, histograms from the selected channels were used to compute the moments of the DTOFs-specifically the sum, mean, and variance moments 79,80 . DTOF moments underwent further preprocessing using a motion correction algorithm, Temporal Derivative Distribution Repair (TDDR) 81 . TDDR can leave spiking artifacts in the presence of baseline shifts-we detected these using moving standard deviation, and corrected them using cubic spline interpolation 82 . Finally, we detrended the data using a moving average with a 100-second kernel.
The relative changes in preprocessed moments were then converted to relative changes in absorption coefficients for each wavelength using the sensitivities of the different moments to changes in absorption coefficients (moments Jacobians, i.e. derivative of each moment with respect to a change in absorption coefficient), following the procedure recently described in 46 . To obtain the sensitivities of the different moments, we used a 2-layer medium with a 12-mm superficial layer. We used a finite element modeling (FEM) forward model from NIRFAST 83,84 , and integrated the moment Jacobians within each layer to compute sensitivities.
Finally, we converted the changes in absorption coefficients at each wavelength to changes in oxyhemoglobin and deoxyhemoglobin concentrations (HbO and HbR, respectively) using the extinction coefficients for the two wavelengths and the modified Beer-Lambert law (mBLL) 85 www.nature.com/scientificreports/ 3. Extraction of physiological signals (a) Pulse rate, pulse rate variability, and breathing rate extracted from TD-fNIRS data: Within each session, we used the sum moment (total counts) of within-module channels (SDS = 10 mm) from the prefrontal plates (prior to performing short channel regression). Data from channels with high scalp-coupling index 86 (SCI > 0.85, a minimum of 5 channels) were averaged. We used the python heartpy module 87,88 to extract physiological signals, for example average and instantaneous pulse rate (PR), and pulse rate variability (PRV). In parallel, we recorded data using an auxiliary PPG sensor and physiological data were also extracted using the heartpy module. The physiological signals computed from the NIRS data were then validated against those measured by the PPG sensor. First, the instantaneous PR from the PPG was interpolated to match the sampling rate of the Flow1-derived PR. For each session, the two PR time courses were compared using a Pearson correlation coefficient. Additionally, the average PR and PRV during a given session were compared against each other and the results are reported in Fig. 2c. (b) Physiological measurements from auxiliary physiological sensors: We applied a fast Fourier transform (FFT) on the time series of CO 2 waveforms measured during the session to obtain power spectra ( Supplementary Fig. 3a). Breathing frequency was obtained by locating the peak of the power spectra during each saline and ketamine dosing sessions ( Supplementary Fig. 3b). SpO 2 and EtCO 2 signals were continuously recorded throughout the dosing sessions and the average value of each time series was computed during each session as additional respiratory features.

Absolute chromophore concentrations
The DTOF is the convolution between the time-resolved temporal point spread function (TPSF) and the Instrument Response Function (IRF). As such, the mean and variance moments of the TPSF can be obtained by subtracting the IRF moments from the DTOF moments. We leveraged calibration data for Flow1 that characterized the IRF for each laser in the system as a function of both temperature and optical power, which allowed us to estimate the IRF throughout the recording session by mapping the instantaneous temperature and power setting to the calibrated IRF value. Absolute moments were then computed by subtracting estimated IRF moments (i.e. mean and variance) from DTOF moments recorded during the experimental sessions.
We then used the analytical expression that converts the mean and variance moments to coefficients of absorption for a semi-infinite homogeneous medium 79,80 . This provides an approximate estimate of absorption coefficients as the semi-infinite assumption is an idealization. These absolute estimates of absorption coefficients can be converted to HbO and HbR concentrations with the mBLL, as described above. To obtain a single time course for HbO and HbR, we computed the median value across well-coupled prefrontal channels with a SDS range of [15,25]  The time series of data from each channel (for each chromophore) was converted to the frequency domain using an FFT. The ratio of the power in the 0.01-0.08 Hz frequency range was calculated relative to the full 0-0.25 Hz frequency range. For time-varying fALFF analysis, we used a window size of 5 min with a stride of 10 s and fALFF was calculated for each window (Figs. 3b, 4a). The fALFF time series was then normalized with respect to the initial value within a given session. For fALFF within the different segments of a session (first half, second half), fALFF was computed over the entire time window of interest (Fig. 3c). The values were then normalized by fALFF computed from a 400 s segment of data in the immediately preceding baseline session (selected to end before the breathing challenge, i.e. from 100 to 500 s) to correct for any differences in baseline fALFF between the saline and the ketamine dosing sessions. Whole-brain fALFF was then calculated as the average over allchannels. (b) RMEQ total score prediction To predict participants' total RMEQ scores, we used three features of PR, HbO fALFF, and HbR fALFF (all time-varying features normalized by their initial value in a given session). We used a Leave One subject Out (LOO) cross-validation combined with linear regression. On each fold, both saline and ketamine sessions of a given participant were left out as the validation set and the RMEQ scores were regressed against the data from all other subjects; predictions on the validation set are shown in Fig. 5b. The quality of the fit was further quantified using the mean absolute error on the validation set (i.e. the absolute value of the difference between average predicted RMEQ and true RMEQ).
To test whether the performance of the model was above chance level, we used a shuffling procedure. Here, the RMEQ scores ("targets") were randomly shuffled across sessions so that the samples for each session were assigned a new RMEQ score from a randomly selected session. The entire LOO CV training procedure was then repeated using these shuffled targets, yielding a coefficient of determination and mean MSE. This shuffling procedure was repeated 1,000 times. We then calculated a p-value for the better-than-chance significance of our models' results by computing the percentage of shuffles which yielded a higher coefficient of determination or lower MSE.
To tease apart the contributions of the neural and physiological features, partial correlations were also computed. Here, both fALFF features were used as covariates to PR, while only PR was used as a covariate for each of the fALFF features. As we are interested in exploring whether the neural signal www.nature.com/scientificreports/ contains additional useful information on top of physiological signals, and since the two fALFF features are highly correlated with each other, including one as a covariate for the other would artificially lower its partial correlation. (c) RMEQ prediction time-series To extract time-varying RMEQ total score predictions, the model (full model, using all features) was evaluated at every time point during the ketamine session. Prediction on the validation set for each participant is shown in Fig. 5c. 6. Functional connectivity (a) Connectivity computation The time series of each channel was further processed by applying a high-pass filter using a bank of discrete cosine basis functions, followed by a low-pass filter using an acausal finite impulse response (FIR) filter. The pass band was 0.01-0.1 Hz, as is customary in the functional connectivity literature. Functional connectivity was computed for all pairs of channels as the Pearson correlation coefficient between their time courses, for each chromophore. For the baseline sessions, we selected a 400 s segment of data ending before the breathing challenge, from 100 to 500 s, to compute connectivity. This segment of data can be considered to be a resting session.
For the dosing sessions, we selected a 400 s segment of data, from 1100 to 1500 s, to compute connectivity. We chose this segment to be near the end of the session (maximal expected neural effect), and of the same duration as the baseline segment. (b) Global brain connectivity (GBC) For each channel in the prefrontal plates, a connectivity strength score was obtained by averaging all positive connections to other channels measured over the whole head. A connectivity strength map (limited to the prefrontal plates) was thus computed for each session.

(c) Statistical Comparisons
We used the baseline session as a baseline for all connectivity analyses. We reasoned that the metric of interest was a change in the state of participants, and using a within-visit reference allowed us to focus on this change. For functional connectivity comparisons, we computed the median of the within-network connections and between-network connections separately for each network and each participant, then used a t-test to assess significance at the group level. For global functional connectivity comparisons, we computed the spatial maps of prefrontal GBC change for the ketamine and the saline sessions for each subject, then used a paired t-test to highlight consistent modulations at the group level. We further looked at the distribution of the GBC changes in the group map (collapsing the spatial dimension) and tested whether the distribution was centered on 0 using one sample t-test. Survey data collection. Participants completed the Columbia Suicide Severity Rating Scale (CSSRS) 89 at the beginning and end of each dosing visit for safety reasons. They also completed the Quick Inventory of Depressive Symptomatology (QIDS-SR) 90 to report depressive symptoms experienced over the previous week at the beginning of each dosing session and during the follow-up phone call. The clinician administered the Brief Psychiatric Rating Scale (BPRS) 91 and Clinician-Administered Dissociative States Scale (CADSS) 92 approximately 60 min after the placebo or ketamine administration. Participants then completed the Revised Mystical Experience Questionnaire (RMEQ-30) 93 and the 5-Dimensional Altered States of Consciousness Rating Scale (5D-ASC) 94 . Participants were also asked about their overall experience with wearing the TD-fNIRS headset using a customized survey ( Supplementary Fig. 1f.).

Data availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.  80-84 (1998).